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We present some exact results for the effect of disorder on the critical properties of an anisotropic 
XY spin chain in a transverse field. The continuum limit of the corresponding fermion model is taken 
and in various cases results in a Dirac equation with a random mass. Exact analytic techniques can 
then be used to evaluate the density of states and the localization length. In the presence of disorder 
the ferromagnetic-paramagnetic or Ising transition of the model is in the same universality class as 
the random transverse field Ising model solved by Fisher using a real space renormalization group 
decimation technique (RSRGDT). If there is only randomness in the anisotropy of the magnetic 
exchange then the anisotropy transition (from a ferromagnet in the x direction to a ferromagnet in 
the y direction) is also in this universality class. However, if there is randomness in the isotropic part 
of the exchange or in the transverse field then in a non-zero transverse field the anisotropy transition 
is destroyed by the disorder. We show that in the Griffiths' phase near the Ising transition that 
the ground state energy has an essential singularity. The results obtained for the dynamical critical 
exponent, typical correlation length, and for the temperature dependence of the specific heat near 
the Ising transition agree with the results of the RSRGDT and numerical work. 
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I. INTRODUCTION 

An important feature of many low-dimensional models of strongly interacting electrons is that they exhibit quantum 
phase transitions, i.e., they undergo a phase transition at zero temperature as some parameter is varied.cl Experimental 
realization of this occurs in heavy fermion materials, such as CeC^-zAu^ which undergo an antiferromagnetic- 
paramagnetic phase transition induced by pressurefllNear the critical point unconventional metallic behaviour is 
observed and is enhanced by the presence of disorder flu It has also been proposed that a quantum critical point plays 
an important role in cuprate superconductors. u Quantum phase transitions in the presence of impurities or disorder 
also occur in 4 He and 3 He absorbed in porous mediaja superconductor- insulator transitions in dirty thin filmsjlthe 
derealization transition in the quantum Hall effect and the metal-insulator transition in doped semiconductors. u 

Compared to thermal phase transitions in disorder-free systems these transitions are poorly understood because 
many of the theoretical methods (e.g., exact solutions, the renormalization group and e expansions) that have proven 
so useful for pure systems at non-zero temperaturestl are difficult to implement for disordered systems.t3 These phase 
transitions are associated with particularly rich physics such as large differences between avecage and typical (i.e., most 
probable) behaviour, new universality classes, logarithmic scaling and "Griffiths phases" ,1111 in which susceptibilities 
diverge even though there are only short-range correlations. Low energy properties of the system are dominated 
by extremely rare configurations oiLthe system. It has recently been proposed that Griffiths phases can lead to 
unconventional metallic behaviour .□113 

Fisher recently made an exhaustive study of the effect of randomness on-what is arguably the simplest model to 
undergo a quantum phase transition: the transverse field Ising spin chainJ13 He used a real space renormalization 
group decimation technique (RSRGDT), originally developed by Dasgupta and MaJiJ which he claimed is exact near 
the critical point. Fisher found the phase diagram (which included a Griffiths phase near the critical point), all of the 
critical exponents (some of which are irrational, as shown in Table |) and scaling functions for the magnetization and 
correlation functions in an external field. It is striking that the latter have never been derived for the disorder-free 
case but can be derived in the presence of disorder because distributions, become extremely broad near the critical 
point. Many of Fisher's results have been confirmed by numerical workJl3~tZI The RSRGDT has now also been used 
to study the effect of disorder on dimerizedllfl ^iid anisotropic spin-i chains, spin-1 chai ns P 9 P °I chains with random 
spin sizes,cil quantum Potts and clock chains,i! and diffusion in a random environment @Ei Possible experimental 
realizations of random spin chains are given in Table |. There is a direct connection between-the critical behaviour 
of the random transverse field Ising chain and random walks in a disordered environment. E3~E3 Fisher's results have 
also been related to the Kondo lattice in one dimcnsion.Efl 
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An important question is whether some of the same interesting physics occurs in higher dimensional models. 
Senthil and Sachdev did find this to be the case in a dilute quantum Ising system near a percolation transition.LJ It 
is particularly interesting that some of the most striking features that Fisher found in the one-dimensional model (a 
variable dynamical critical exponent which diverges at the critical point and the average and typical correlations are 
associated with different critical exponents) have recently been found in the two-dimensional random transverse field 
Ising model.EJ 

The outline of the paper is as follows. In Section II we introduce the model, an anisotropic XY spin chain in a 
transverse field, where all the exchange integrals and transverse field-are random. A similar model was also recently 
studied numerically using the density-matrix renormalization groupH A Jordan- Wigner transformation is then used 
to map the model onto a non-interacting fermion model. Section III contains a brief summary of the known properties 
of the disorder- free model that are needed to understand the rest of the paper. In Section IV we take the continuum 
limit of the fermion model for various cases. The Ising transition and the anisotropy transition with only randomness 
in the anisotropy that results in a Dirac equation with a random mass. The isotropic XX chain in a transverse 
field with randomness in the exchange and/or transverse field reduces to a Dirac equation with a random complex 
mass. Mapping the spin chain to these Dirac equations has the advantage that a number of different exact analytic 
techniques can then be used to evaluate the density of states and the localization length. The properties of the solutions 
corresponding to the universality class of the random transverse field Ising model are then discussed in Section V. 
By examining the energy dependence of the density of states we evaluate the dynamical critical exponent, show the 
existence of a Griffiths' phase near the transition, and show that the ground state energy has an essential singularity 
at the transition. We also present results for thermodynamic properties and the typical correlation length. The 
results obtained for the dynamical critical exponent, typical correlation length, and for the temperature dependence 
of the specific heat near the Ising transition agree with the results of the RSRGDT and numerical work. Since our 
approach is explicitly exact our results are consistent with Fisher's claim that the RSRGDT gives exact results for 
critical properties. In Section VI we point out that the properties of the incommensurate solution are such that it 
implies that the anisotropy transition is destroyed in a non-zero transverse field if there is randomness in the isotropic 
exchange or in the transverse field. A brief report of some of the results presented here appeared previouslyB3 



II. THE MODEL 



The Hamiltonian to be considered is that of an anisotropic XY spin chain in a transverse field: 

L 

H = - E ( J n«+1 + Jl°K+l + K<) • (1) 

The a", (a = x, y, z), are Pauli spin matrices. This is a quantum model because the Pauli matrices do not commute 
with one another. The interactions, J x , J^, and transverse fields, h n , are independent random variables with Gaussian 
distributions. All the results given in this paper are for this ferromagnetic case but also hold for the antiferromagnetic 
case. By means of a spin rotation J x and h n can always be chosen to be non negative. We shall assume that J% is 
also non negative so that there is no frustration in the system. The average values will be denoted 

{J x n )^J\ {Jl)^J y , {K)=h. (2) 

The deviation of the random variables from their average values is assumed to be small, relative to the average value. 
We can write our three parameters in terms of a random part and an averaged part, 

J a n = J a + 5J% (a =x,y), h n = h + 5h n . (3) 

The random variables are uncorrelated between sites and have variances 

<(J°- J a ) 2 ) = (SJ a ) 2 (a = x,y), (( hn -h) 2 ) = Sh 2 . (4) 

For jy = the model is the random transverse field Ising spin chain which is the quantum analog of the two- 
dimensional Ising model with random coupling in one direction, introduced by McCoy and Wu,E3 and studied by 
Shankar and Murthy.Efl 

At zero temperature and in the absence of disorder the model undergoes two distinct quantum phase transitions .EEI 
Both transitions are second order. The phase diagram is shown in Fig. [lj The transition at J x + J v = h from a 
paramagnetic to a ferromagnetic phase will be referred to as the Ising transitions^ The transition at J x — J v for 
h < (J x + J y ) from a ferromagnet with magnetisation in the x direction to one with magnetization in the y direction 
will be referred to as the anisotropic transition£Ix3 This paper considers the effect of disorder on these transitions. 
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A. Mapping to a fermion model 

We perform a Jordan- Wigner transformation which maps the Pauli spin matrices in ([!]) onto spinless fermionsl 3 ^' :l 
The Pauli spin matrices 



01\ v ( -i\ z I 
10 = ( i ) * = -1 



(5) 



satisfy the algebra 

la a ,a b } = 2ie abc a c (a a ) 2 = 1. (6) 
Define the following new operators on each site, 

4 = \{< + i<) «" = \K - to*). (7) 

The inverse transformation is 

<rn = a i + a n er* = i(a n - a* ) <?n = 1 ~ 2a n at. (8) 

Using these definitions the following relations can be obtained, 

{a\, a n } = 1, a 2 n = a\ 2 = 0, [a f a„] = [a} n , ajj = [a m , a„] = m^n (9) 

which show that the a„'s and ajj's are neither fermion operators nor boson operators. The Hamiltonian, (Q), in terms 
of these new operators is, 

L 

H=-J2[(Jn + J v n){aia n+1 + ai al +1 ) + (J* - J^)(a+4 +1 + a n a n+1 ) + h n {l - 2a n aJ, )]. (10) 

n=l 

Now consider a second transformation, 



c„ = exp(7ri a]aj)c 



j'=i 

n-l 

4 = 4 cxp(-vri ^ a] aj ). (11) 
i=i 

The c n 's and c]j's are fermion operators satisfying the following anti-commutation relations, 

These fermions can be viewed as kinks or domain walls in the local magnetization.0 The Hamiltonian is now, 

L 

H = -J2^ J n + J n)( C i C n+l ~ C n C ] n+1 ) + {J* - J^)(c„4 +1 - C n C n+1 ) + h n {c\ lCn - C„C„)]. (13) 
71=1 

The boundary terms have been neglected since they do not contribute to the thermodynamic limit. 

III. SOLUTION OF THE DISORDER-FREE CASE 



The model in the absence of disorder has been solved previously.HciEj We now highlight certain aspects of the 
solution that will turn out to be particularly relevant to the effect of disorder. In the disorder free case J* = J x , 
Jf % = J v and h n — h so that 

L 

H = -J2i( JX + J v ){cic n+1 - c n cl +1 ) + (J x - J y ){cicl +1 ~ c„c„+i) + h(cic n - c„ C t )]. (14) 

n=l 
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The case h = corresponds to the anisotropic XY spin chain and was first solved by Lieb, Schultz and Mattis.tl The 
case J v = is the transverse field Ising chain and was first solved by Pfeuty.Ej 
We introduce the Fourier transform of the fermion operators 

V L k 

c ^ = ^E c Ie- mfc - (15) 
Periodic boundary conditions (c„ = c„ + i) require the wave vector, fc, to take the following discrete values, 

* = ^p> m = ~L,... ,0,1,. ..,^L-1, (16) 



assuming L to be even. Substituting (|15|) in the Hamiltonian ([14]) gives 

H = - £)[2(( JX + JV ) cos k + h ) c W +i(J x ~ J v ) sin fe(4 c -fc + c kC-k) - h] . (17) 

k 

This Hamiltonian maybe diagonalised by the Bogoluibov transformation, 

cj, = cos <t>{k)b\ + i sin cf)(k)b-k 
Ck = cos 4>(k)bk — i sin cf)(k)b { _ k , (18) 

where the b^s and b\js are operators with fermion statistics, and 

,„,,,,, (J x — J v ) sin k , „, 

tan(2^(fc)) = /T ; , TfA ' , ■ r - (19) 



This gives, 



where 



(J x + Jy)cosk + h 

H = -Y,E(k)[blb k - 1/2] (20) 
fe 

£(fc) = 2(/i 2 + (J 1 - J^) 2 + 2/i( J* + J») cos k + AJ x J y cos 2 fc) 1 / 2 , (21) 



and we have used, cos fc = 0. 

The energy gap will occur at a wave vector kg such that 

mk^ = Q 

afc 

We shall always take fco to be the positive solution of the above equation. Because the energy is symmetric in k there 
will be two energy gaps at ±feo- To describe the solutions of this equation it is convenient to define 



h(J x + jv) 
a = , 

u x jy ' 

which is always negative for non negative J x , J v and h. The energy gap wave vector, ko, is given by 



(23) 



cos kg = a a > — 1 , 

k = n a < -1. (24) 



Typical dispersion curves for these two cases are shown in Fig. 0. If we define, 



_ jy 

1=- — , (25) 

j x + jy y ' 



and express a as 
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then, the boundary between the two cases may be defined by 

h 



J x + jy 



= 1 - 7 ■ (27) 



This boundary is shown as a dashed line in Fig. |l|. The two cases correspond to a commensurate (a < —1) and an 
incommensurate (a > —1) phase.Ea It will turn out that the effect of disorder on these two phases is very different. 
The energy gap at k — fco is 2 A where 

A = E(ko). (28) 

The system is at criticality when the gap vanishes, A = 0. When a < — 1 (fco = 7r) the gap vanishes along the line 
h = J x + J v . We shall refer to the corresponding phase transition as the Ising transition. When a > — 1 the energy 
gap vanishes along the line J x — J v , providing h < J x + J v . We shall refer to this transition as the anisotropic 
transition. The lines along which the gap vanishes are shown as solid lines in Fig. ^. 

The critical behaviour is determined by those low energy states near the energy gap where k ~ fco. If fc — fco is small 
the energy can be written as a Taylor series, 

E{k) 2 = A 2 +v 2 (k-k a ) 2 + ... (29) 

where 

v a = 2[AJ X jv - h(J x + J y ) cos fc - 8 J x J y cos 2 fc ] 1/2 

A = 2[h 2 + (J x - J y ) 2 + 2h(J x + J y ) cos fc + 4J x J y cos 2 fco] 1 / 2 . (30) 



A. The ground state energy 



The ground state energy, e, of the Hamiltonian is the energy of the filled Fermi sea 

r * dk 



£ (A) = - 



f dk .,. 



(31) 



where E(k) is given by (J2l|). To discover the nature of the singularity at A = we differentiate the above integral 
with respect to A 2 . For small A the differentiated integral is dominated by those low energies states close to the 
energy gap. Hence, we need only consider those states determined by the low energy dispersion relation (E9), 



de 



fco + fec 
kg— k c 

1 

27TUo 



dk 



1 



4^ ^A 2 + Uo 2 (fc-fco) 2 
A 



In 



2w fc c 



where fc c is a cut-off wave vector. Integrating with respect to A 2 gives 



c(A)-e(0) 



A 2 



l-21n 



A 



2w fc c 



(32) 
(33) 

(34) 



The singularity of the ground state energy is thus logarithmic. The critical exponent a, defined by e(A) ~ A 2_Q , is 
a = + . This critical exponent corresponds to the specific heat critical exponent of the corresponding two-dimensional 
classical Ising model. 
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B. The magnetization and correlation length 

Barouch and McCoyEi!l calculated the magnetization and correlation functions Jar the disorder-free model. Its et 
al. considered the case 5 = 0.c3 Further analysis was done by Damle and Sachdev.cil The magnetization, M x = (a x ), 
and the correlation length, £, are defined by the asymptotic behavior (r — ► oo) of the correlation function, 

(« l+r > -> {M x f + £ exp(-r/0 (35) 

where A is a constant. If h > J x + J v the system is a paramagnet and the magnetization is zero. If h < J x + J v and 
J x > J v the system is a ferromagnet in the x direction, and the magnetization is, 



rx\2 , nr 2 7' 



(M X Y = (-l) r 



1 + 7 



1 



J x + jy 



(36) 



This implies that the critical exponent (3 is 1/8 for the Ising transition (approaching the transition as a ferromagnet) 
and 1 /4 for the anisotropic transition. 
The correlation length, £, is given by 

exp(^)HA 2 r 2 (37) 

where, 

h ~|2 _ ft _ „,2^i 



[b-£nr) a -(l-7 



A 2 = " (38) 

This quantity is real (complex) in the commensurate (incommensurate) phase. As a result 

f— 1\ f A^ 2 commensurate 
6X ^ \ £ y \ -jq^ incommensurate. ^ ' 

This implies that the critical exponent v — 1 for both the Ising and anisotropy transitions. 
On the Ising critical line 

(«+r) - ^4 ( 40 ) 

and the critical exponent r\ = 5/4. On the anisotropic critical line 

(«+r) - ^Jj2 ( 41 ) 

and the critical exponent r\ = 3/2. The anisotropic transition has the same critical behavior as a pair of decoupled 
Ising models. The critical exponents for the Ising transition are summarized in Table |l[ 



IV. THE CONTINUUM LIMIT 



We shall now look at the effect of disorder on the critical behaviour. To do this we take the continuum limit of 
the disordered Hamiltonian written in terms of Fermi operators, (|l3|), when the system is near criticality. We will 
assume that, for weak disorder, the phase transitions of the disordered system are close to the phase transitions of the 
disorder free system. Those fermion states most effected by the addition of disorder will be those low energy states 
near the energy gap, that is, those with wave vectors near ±/co- 

The Hamiltonian may be broken into slowly and rapidly varying parts. This is done by replacing the Fermi operator, 
c„, with two slowly varying functions, ijjR(n) and - 0i( ri )i which describe right and left movers respectively, 

c n = ±=[e-^ n 4> R {n) + e tk ° n ^ L ( n )}. (42) 

The exponential terms represent the rapidly varying part of c„. From the anticommutation relations of the Fermi 
operators, c„ and cj,, it is possible to derive anticommutation relations for the slowly varying functions, 

W R {n),^R{m)} = {ijjl(n),ip L (m)} = 5 nm , (43) 
and all other combinations are zero. 
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A. Ising Transition 



When k — ±7r, equation ( |42| ) may be simplified to, 

c„ = (44) 

with, 

{V t (n),V(m)} = 5„ TO . (45) 

After substituting this into the Hamiltonian, we take the continuum limit. To do this we transform the discrete 
variable, n, into a continuous variable, x, and we write, 

ip(n) = ijj{x = n), ip(n + 1) = ip(x) + d x ip(x). (46) 

The function ip is slowly varying and so its derivative is very small. Where appropriate we can neglect these derivative 
terms. The disorder is assumed to be small, hence terms combining both derivative terms and disordered terms may 
be neglected. Substituting these approximations into the Hamiltonian and replacing the sum over n with an integral 
over x gives, 

H = ^2(tf,il>)[iJ-o»d a + (J + (x) - h(x))a*} (t), (47) 

71=1 ^ ' 

where 

J±{x) = J x n=x ± J% =x , h(x) = h n=x . (48) 
By performing the following rotation into a new set of Pauli spin matrices, with ko — ±7r, 

AJ + sinfco 1 vo(h + J + cos k ) 



-a* + - — ~r~~. ; — t(7 3 



2J-(J+ + hcosk ) 2J-(J+ + hcosko) 

A cos kg 1 vq sin ko 3 
rC + „ , T ; ; — r(T 



2(J+ + h cos k ) 2( J+ + h cos k ) 
a x =a 2 , (49) 



where A and vo are defined equation (p0|), the Hamiltonian becomes, 

H =\ J dx^(x)^[-iv a 3 d x + V{x)a 2 ]^{x), (50) 

where, 

V(x) = 2\J+ -h\± 2{6J+(x) - Sh(x)), (51) 

and SJ + (x) and 5h{x) are the random parts and J + and h are the average parts of J + {x) and h(x) respectively. The 
function V(x) is real and,its average value is A, the energy gap of the pure system. For the case of no disorder (|50| ) 
was derived by Shaipkar.tJ The case of the transverse field Ising chain with randomness only in J x or h was derived by 
Balents and Fisher O The fact that the Ising transition is described by the same equation for any anisotropy shows 
that it will be in the same universality class as the random transverse field Ising chain studied by Fisher. 



B. Anisotropic Transition 



Near the anisotropic transition we must use the more general decomposition of the Fermi operators shown in 
equation ([42|). As was done near the Ising transition we replace the discrete variable, n, with a continuous variable, 
x, and replace discrete differences with derivatives. Next, we remove those disordered terms which are negligible. We 
shall make one approximation which was not necessary near the Ising transition. We will neglect all rapidly varying 
terms. A rapidly varying term may be neglected because its integral will vanish. Terms involving the product of two 
rapidly terms may not be neglected since the two rapid variations may produce a slowly varying part. The exponential 
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terms are rapidly varying and the random terms may have slowly and rapidly varying parts. Combining all these 
approximations gives the following Hamiltonian, 



H = - dx 



i(J + (x) cosfco + h{x)){i^ R ijj R + ip^ipL) + iJ (x) sinfc o -0i-0lj - 



+e 2tk ° x (6J + {x)e lk ° + 5h(x))il>y; L ] + h.c. (52) 



When kn = —j, J x (x) — J y {x) — and h(x) — 0, this Hamiltonian is equivalent to that obtained by Balents and 
Fisheries The Hamiltonian in equation ( p2[ ) can be simplified in two particular cases. 

1. Case I (6 J* = -5J%, 5h n = Q) 

If the only randomness is in the anisotropy, which must be equal but opposite (5J + (x) — Sh(x) — 0), the Hamilto- 
nian reduces to, 



H = - J dx{ij)^ R , iPl)[(J + cos/cq + h)a z - J (x)sink <r y 



+i{J+ sin k (J z + J- cos k <T y )d x ] ( V (53) 
On performing the rotation in equation (ff9|), the Hamiltonian can be simplified to, 

H = \J dx(i> R ,i> L )[-iv * 3 d x + Vix)* 1 } (tf\ (54) 

We have used the definition of fco on the anisotropic critical line, cosfco = — -jt, and defined, 

V(x) = 2\J-\smk ±2SJ-(x)smk . (55) 

Note that, like the Ising transition, V is real and its average value is A. Hence, in this case the anisotropy transition 
is in the same universality class as the random transverse field Ising chain. 

2. Case II (J% = J%) 

Another special case of the Hamiltonian in equation ( |52] ) is when the is no anisotropy (J~(x) = 0). Note that, 
in the absence of disorder this restricts the model to the anisotropic critical line. It shall also be assumed that any 
disorder is rapidly varying. The new Hamiltonian is, 

H = -\ J{i>i,i>{)[-2iJ + sink Q cx z d x + i{xYc7 + +i{x)a-] , (56) 

where = 2e~ 2lk ° x (SJ + (x)e~ lko +5h(x)) and = \{o x ±ia v ). Since fco is incommensurate with the lattice £(x) 
is complex. We have neglected a term involving the sum over the magnetic field since it is a constant. Consider the 
following rotation, 

d = —a 
a = a 

a 2 = -a y . (57) 

With this rotation, 

H = \f *{xy[-iv a*d x + Z{x)<j + +£(aOV-]¥(aO, (58) 
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where a ± = ^(er 1 ± id 1 ). The complex function, £(x) = 2e 2lk ° x (d J (x) + e lk ° + Sh(xj), has the following properties, 

(£(*)) =0 

(e(^'))=o 

(^)^'r)=7^-^) 5 (59) 

where 7 = 4((<5J+) 2 + (^) 2 ). 

In summary, all three Hamiltonians, (|50|), (Q) and (|5§|), can be written in the form,c3 

H =\ J dx^{x)^ [~ivo<J 3 d x + V{x)a+ + V(x)*a~]^(x). (60) 

Note that the structure of the spinor, 'J, is quite different in all three cases. The function V(x) satisfies, 

(V(x)) = A 

(V(x)V(x')*) = A 2 + jS(x — x 1 ). (61) 

The anisotropic case II ( J~ = 0) has A = 0. The fact that V{x) is complex for the anisotropic case II will lead to 
qualitively different behaviour. In fact, in that case the disorder removes the phase transition. We refer to the case 
where V(x) is real as the commensurate case (that is, the Ising transition and the anisotropic case I). The casa-where 
V(x) is complepc-is the incommensurate case. The case of real V(x) also describes dimerized XX spin chainsEJcJ and 
spin ladders.^ The case of complex V(x) also describes an XX spin chain in a transverse field with a modulation 
of the exchange with wavevector 2/co- 

V. EXACT SOLUTIONS 

It is useful to define an energy D and a dimensionless parameter 5 which are measures of the disorder strength and 
the deviation from criticality, respectively 

7 A 

D=-L S ^ 6 2 
v D 

Note that for the Ising transition with, J y = 0, to leading ordex-in A/J x , for a Gaussian distribution this parameter 
8 agrees with the 5 defined by Fishei£3 and Young and RiegerJIS 

<ln/*}-(lnJ*> 



5q,pi 



~ ((hi ft) 2 ) - (Inh) 2 + ((In J x ) 2 ) - (In J*) 2 ' 

The advantage of casting the problem in the form of the Hamiltonian (^) is that the latter has been studied 
extensively previously, and exact analytic expressions given for the energy dependence of the disorder-averaged det 
of states (p(E)) -anxLlhe localization length^ E) . The exact results-have been found by Fokker-Planck-equationsJ: 
supersymmetry,E3oo the replica trickJ£3~E3 S-matrix summationJ£3 and the Dyson-Schmidt method.Ej 

Due to the one-dimensionality all the states are localized by the disorder. The localization length can be found 
because in one dimension it is related to the real part of the one-fermion Green's function.EjCj 

The density of states and the localization length are related to the one-electron Green's function G(x, x, E) and 
can be written in terms of fg(u), the derivative of a dimensionless function fg(u), 

Tr G(E,x,x) = A J_ +in P (E) = wp f' s (E/D) (64) 

where po = 1/(ttvq) is the value of the density of states at high energies (\E\ ^> A,£>). The function fg is different 
for the commensurate and incommensurate cases, 



i — u-S- \ii\H x (u)] commensurate /r , r \ 

/*(") = i xdTn «V • + ( 65 ) 

I o-gg m[li u {d)\ incommensurate 

(2) 

where Hg is a Hankel function of order 5 and Ii y is a modified Bessel function with imaginary index. 
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A. Solution using a Fokker-Planck Equation 



To demonstrate how an exact solution may be found we will derive the density of states for the commensuraite[eap« 
(V(x) real) by using Fokker-Planck equations. Many authors have studied mathematically equivalent systems .E3c3 : E2l 
Consider a general Dirac type equation, 

dtp 

- iv 0-g^ + V(x)l/)2 = Elpl 

iv ^ + V(x)ip 1 =Eip 2 . (66) 

The function, V, is a real and random function of the form, 

V(x)=A + £(x), (67) 
where A is a constant and £ is a random field which obeys the following statistical averages, 

(£(*)> = 0, (Z(x)Z(y)) = >yS(x - y). (68) 
We reduce the Dirac equation into a system of equations for two real functions by the following transformations ,0 



* \ ( vi + r 2 \ ( * A _ (a - r 2 



(69) 



then we let, 

It can be shown that (/i, /2) and (0i, ^2) satisfy the same equations, 



v £ -V(x) J\hJ~ \f: 
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Define the following function, 



fx 



E[ J l). (71) 



(72) 



By differentiating the dynamic variable, z, with respect to x and using equation ( |71|) a dynamic equation for z may 
be constructed, 

v ^- = -(E-A)-z 2 (E + A)-az 2 -l)- (73) 
This equation allows us to write down a Fokker-Planck equation for the random variable, z, 

The function P(z,x) is the probability density distribution function of the random variable z at the point x. It is 
the derivative (with respect to x) of the probability that z is less than x. Since the probability of z being less than 
infinity is unity we expect, 

poo 

P(z,x)dx=l. (75) 

This is an important concept when dealing with probability densities. 

We create a stationary Fokker-Planck equation by taking the limit as x goes to infinity. The limit of the probability 
density is, 

lim P(z,x)—p(z). (76) 
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The stationary Fokker-Planck equation is 

1 d 



= 



vq dz 



(E - A) + z\E + A) + j(z 2 - l)^(z 2 - 1 i ) pi-)- 



(77) 



This equation can be integrated and Ovchinnikov and ErikhmanEl showed that the constant of integration is N(E), 
the number of states below the energy E, 



1 



''0 



D 



d 



N(E) = -- ( (E - A) + z\E + A) + - 1)^^ - 1) ) p(z) 



(78) 



To simplify the solution of the above differential equation we perform two transformations. 
The first transformation is defined by the following function, 



cot 



a(x) 



fi 



(79) 



Because of the nature of the cotangent, we can obtain all possible values of the ratio, — /2//1, by restricting a to some 
interval of length 2tt. We will restrict a to the interval [— ■§, 4f ]■ We are only interested in the case where x — + 00. 
As in equation (f75|), the integral of all probability densities of a at large x, p(ct), must equal unity, 

p(a)da = 1. 



The relationship between p(z) and p(a) is, 



The second transformation is, 



p{z) = —2p(a) sin' 



cos a = ± sech </ 
sin a = ± tanh < 



a 



(80) 



(81) 



(82) 



The upper sign refers to a 6 [— -| , -| ] and the lower sign, a £ [|, ^]. It can be shown that p(a) = p(<p) cosh^ so that 
equation (|o|) becomes, 



P + {4>)d<j ) + / p_(0)d0 = l, 



(83) 



where the ± subscript on p refers to the different signs in equation (|S2|). After performing the two transformations 
equation (^sh becomes, 



2£; cosh 6 ± 2 A + 2D 



d_ 

80 



p(<f>)=v N(E), 



This first order differential equation can be solved with the boundary condition that p((f>) vanishes as cf> — > 00, 

v N(E) 



2D 



dx exp 



E A 
— (sinh (j> — sinh x) ± — ( 



(84) 



(85) 



To find the number of states we recall that the probability density must be normalized so that the integral over all 
possible values of the random variable <f> is unity. After some rearranging, 



/oo />oo 
d<t> j da; exp 

dx exp 



E A 

— (sinh <f) — sinh x) + — ((f) — x) 



(sinh i 



sinh a;) — —(<j) — x) 



(86) 



If these two integrals are combined and we let 2y = x — <f>, S = A/D and u = E/D we obtain, after changing the 
order of integration, 
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2Dv~ 1 N(Ey 1 = 4 / dy cosh2Sy #exp[tt(sinh</> - sinh(2y + <j>)]. (87) 

JO J-oo 

Now we let z = <f> + y, 

/>oo />oo 

2Dvq 1 N{E)- 1 = 4 / dycosh(25y) dz exp[— 2m sinh y cosh z] 

Jo J-oo 

/>OC />OC 

= 8 / e?y cosh(2<5y) / dzexp[— 2usinhy coshz] 
Jo Jo 

/•OO 

— S dy cosh(2Sy)Ko(2usmhy) 
Jo 



= 7T 2 [J s (u) 2 +Y 5 (u) 2 } 

where Js(u) is a Bessel function of index S and Ys(u) is a Bessel function of the second kind of order S. The number 
of states with energy less than E is, 



2D 

To find the density of states, p(E), we differentiate N(E) 



N ^ = -,.. It ,ev mV ,e w - (89) 



n(F , 4 J s (u)Ji(u)+Y s (u)Y>(u) 

P^l = 2 \T I \2 > V ( \212 ' ^) 



This is the density of states from which we obtain our results. 



VI. PROPERTIES OF THE COMMENSURATE SOLUTION 



A. The density of states 

Fig. [| shows the energy dependence of the density of states from equation ( pp| ) for a range of values of 5 with 
po = l/(nvo). In the low energy limit (small u) we can take the following approximation of equation (|9(i|), 

£ffi = _12W (91) 

PO 7T Ys[u) 6 

since the Bessel functions of the first kind remain finite for small u whereas those of the second kind become infinite, 
as shown in the following small u approximations, 

W~{:i r 2 (j)(i „ r < (92) 

These low energy limits are substituted into equation (^Tj) and then the dominant terms are retained, that is, the 
smallest powers of u, 



P{E) f B[ln(*)] 



2D " 2S-i P , „ (93) 



The divergence in the density of states at E — is sometimes referred to as the Dyson singularity. 

The function Ys(u) is continuous as u and 5 approach zero. This property is not apparent from equation (93). To 
avoid this problem we take the small 5 (close to criticality) limit of equation (|9l]) before we take the small u limit. 
Before we take any limits we write the Bessel function of the second kind in terms of the Bessel function of the first 
kind, 

J s (u)cq S (6k)-J- S (u)_ 
sin (on) 
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We can find a small 5 approximation to J$(u) from its series expansion, 



J s {u) 



S oc 

£ 

fe=0 



h 2 ) 



kW{5 + k+ 1) 



u) J {u). 



(95) 



So now we have a small 6 approximation for a Bessel function of the second kind, 



Y s (u) = J (u 



auY-auy 



Sty 



When u is small we set Jq(u) = 1 and equation (M) becomes 



^ = 2^ (— 



2(5-1 



i-(^) 



2-> 



(96) 



(97) 



This agrees with the scaling form obtained by Balents and Fisher.0 By taking appropriate limits it can be shown 
that this formula agrees with equation (p3[). For 5 — we use, 



Urn 

<5^0 



E \ 25 
2D) . 



lim — 

8^0 1 



1 



(1 



2S\n^] 



and for 5 becoming small we use, in equation (|9l|), T(S) ~ l/<5. 

The low-energy (|S| <C D) dependence of the density of states contains some important physics. The density 
of states diverges at E = for 8 < 1/2 and is zero at E = for <5 > 1/2. These two cases lead to qualitatively 
very different behaviour. In the former case some susceptibilities will diverge as the temperature approaches zero. 
This corresponds to a Griffiths or weakly-ordered phaseO Hence, for the Ising transition there will be four phases: 
ferromagnet, weakly-ordered ferromagnet, weakly-ordered paramagnet, and paramagnettj (see Fig. ||). 



B. Ground state energy 



The dependence of the ground state energy of the disordered commensurate system on 6 can be shown to be 
infinitely differentiable, but nO|t-analytic. To show this we follow a procedure similar to that used by McCoy and 
WuE3 and Shankar and MurthyEil who considered the analogous two dimensional classical system. To find the ground 
state energy in the presence of disorder we use, 



e(S) 



p{E)EdE. 



(99) 



We make use of the expression (|97]), for the density of states at low energies, which we assume is accurate up to an 
energy E c , which is less than 2D, 



e(S) = -np D 



2 Q r3 



Z?2(5 



E 26 (l + E 2d ) 
(l-£ 25 ) 3 



dE 



p(E)EdE. 



By integrating by parts, 



8<f 



3 r2n E 2 \l + E 



2ti\ 



(1 - E 2S f 



-dE 



45 2 E 1+2S 



[l-E 2S f 



2SE 



(1 - E 2S ) 



25 



dE 



l_ E 2&- 



(100) 



(101) 



As for the disorder free case (see Section III A| ) we subtract off the ground state energy at 5 = 0. To calculate the 
S = ground state energy we require the limit in equation (|98|), with which we obtain, 



E r 



E r 



e(0) 

irpoD 2 2D(\n^) 2 ^ 2D\n^ 



dE 



lim 



1 



InE S^o TTp D 2 J E 



p(E)EdE. 



(102) 
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By subtracting the zero 5 case from the small 5 case and combining those terms analytic in 5 in a function /(£), we 
obtain, 



e(S) - e(0) = -Trpo-D 2 



-npoD* 



dE 



25 



1 



1 -E 2S 
d^ 2S [(l- 



bxE 



f(S) 



(103) 



where we have made the substitution E = e £/ 25 and set E c = 2D, because-this does not affect the analytic properties 
of the integral. This integral can be solved in terms of Euler's tp function,Eil 



1 



e(S) ~ e(0) = TT Pa D z In 25 + 7rp D^ I — J + f{5) 

OO 

= -irp Q D 2 5 - irp D 2 £ B 2n (25) 2n (2n)- 1 + f(5), 



using the small <5 approximation for the ip function.El The Bernoulli numbers B2 n are proportional to ^"L for large 
n. Because of this, the ground state energy has zero radius of convergence about the point 5 = 0. Thus the ground 
state energy is infinitely differentiable but is not an analytic function of 5. The critical exponent, a, defined in Table 
U cannot be defined in this case. 



(104) 



C. Thermodynamic properties 

1. Free energy 

For any particular configuration of the disorder the free energy per site of the system is 

F = -fc B T^lnf2coshf^)) (105) 



where {Ek} denotes the eigenvalues of the Hamiltonian (|60|). This simple formula holds because the eigenstates of 
the Hamiltonian are non- interacting fermions. It then follows that the disorder- averaged free energy is 



(F) = -~k B T dE(p(E)) In ^2 cosh 



E 



2k B T 



(106) 



The low-temperature behavior of the specific heat and the transverse susceptibility (for the anisotropic transitions) 
follows from the energy dependence of the disorder-averaged density of states£3'E2l We now show this in detail. 



2. Specific heat 

The disorder averaged specific heat is, 

(C( T )) = -r ^P = 1 jf° dEE 2 (p(E))^L (107) 

where f(T) is the Fermi distribution function. For the commensurate case with 5 the mean specific heat is, 
, r(T ., 2n5 P0 f E E 28 + l e -E/ kB T x .oo Q f 

{C{T)) = fc^TlPP Jo dE (l + e-Wy + t ] Ec dE P^QE 

2-n5o n k 2S+1 T 2S r Ec/kBT y 2S+i ~ y \ r°° Bf 

1 dy l ■ * +1 / dEp{E)%. (108) 



T(5) 2 (2D) 2S - 1 J Q y (l + e-y) 2 T J Ea ry 1 dE 
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As the temperature becomes very small the limit of the integral, Ec/ksT, becomes very large- The first integral will 
dominate the specific heat and will be evaluated from zero to infinity. By using integral tablesEH it can be shown that, 

y _w_ , ^v- l > x> _ x 



r°° v x e~ v v-^ 
/ d VT ^— — ^=r(a;+l)V 



k x ' 

fe=l 

= T{x + \){\-2 l - x )C{x). (109) 
where £ is the Riemann zeta function. The mean specific heat is then, 

{C(T)) = 4n6pok B DT(26 + 2)(1 - 2~ M )C(2J + 1) /fcgT\ 26 _ (n0) 

For small <5, (C(T)) ~ <5 3 r 215 , in agreement with FisherEl and the numerical work of Young.0 
The specific heat in the commensurate case with S = is, 

iriTW nD P° rITP Ee- E '^ T 1 9/ 

<C(T)> = -k^L ^ [ln(^)]3(l + e -^T )2 T 7 Ec 

E a /k B T _« 

ye y 



[l n (^)]3 (1 + e - y) 2 



(111) 



To simplify this equation we note that the term ye v /(l + e y ) 2 is appreciable only for values of y of order unity. 
Since we are taking a low temperature limit, T <C D, we may approximate the term, (In vk ^ T ) 3 , to simply (In ^jy-) 3 
when y ~ 1. As in the previous case, the second integral is negligible as the temperature approaches zero and the 
limits of the first integral are zero and infinity. The specific heat for 5 — is 

V V " (m^f) 3 



which has the same temperature dependence as found by Fisher 



3. Transverse susceptibility 

The mean transverse field susceptibility is, for the anisotropic transition 

<X«Cn> = jH dE(p(E))^L. (113) 

For the commensurate case with S ^ and the temperature approaching zero the calculation of the susceptibility is 
similar to the calculation of the specific heat with 5 =/= 0. The mean transverse susceptibility is, 

^^ ^wwM-ff-D ^"- (114) 

If 6 < i and T — > the susceptibility becomes infinite. This is the Griffiths phase r egion . We cannot take the limit 



as 5 goes to zero of the susceptibility since the condition on the integral in equation (109) is x > — 1 which means, in 
this case, S > 0. If we try to take this limit we see that it does not exist. The critical exponent, 7, is not defined. 

Similarly to the specific heat calculation with S = it can be shown that when 6 = the mean susceptibility at 
low temperatures is, 

(Xzz(T)) = / dy i^rT^r; -• (115) 

The integrand is large at y = 1 and y = so we can approximate the integral to an integral from zero to some cut 
off, A, of order unity. We also notice that, for y <C 1, 
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y( l n [^])3 (1 + e -, )2 „(ln[**sT])3- 
Using these approximations the susceptibility is, 



(116) 



(X«(T)> = - 



1 



2/(ln[^]) 3 



k B T[H^T)f 



(117) 



This susceptibility is finite when 5 — unless T = 0. A comparison with the 5^0 result shows that the susceptibility 
is not continuous at the phase transition which is at S = 0. 



D. Dynamic critical exponent z 

This relates the scaling of energy (or time) scales to length scales. We can make the following crude scaling argument 
to extract z from the low-energy behavior of the density of states. The total number of states (per unit length) with 
energy less than E, N(E) scales with the inverse of any length scale £. By definition E ~ l~ z . This implies that 
(p(E)) ~ E x / Z ~ x . Thus for the commensurate case, to leading order in 5, 

Z = W (118) 

in agreement with the renormalization group results of FisherEl and the numerical results of Young and RiegerEl 
This is a particularly striking result because it shows that (i) z ia.not universal and (ii) z diverges at the critical point. 
The latter implies logarithmic scaling and activated dynamics. o 



E. Finite size scaling 

Monthus et aZ.EH studied an equation equivalent to j6C| ) with V{x) real and A = 0.H They have shown that on 
a line of length L, for a typical potential V(x) the lowest eigenvalue E scales like Eq ~ exp(— cL 1 / 2 !^ where c is a 
constant. This is consistent with the scaling pf \uEq with L 1 / 2 at the critical point found numerically S3 The average 
(Eq) ~ expLr-dL 1 ' 3 ) where d is a constant,t£l showing the discrepancy between average and typical values. Fisher 
and Youngc3 recently derived the distribution function for the energy gap from the RSRGDT and compared it to 
numerical results. The distribution function they derived gives average and typical values in agreement with the above 
results. 



F. Correlation lengths 

FisherE! stressed the distinction between average and typical correlations. If = (AiAj) denotes a correlation 

function of a variable Ai then the average correlation function C av (r) = j- Yli=i Ci,i+r is what is measured experi- 
mentally. Away from the critical point C av (r) ~ exp(— r/£ av ) where £ a „ is the average correlation length. However, 
C av (r) is dominated by rare pairs of spins with Cy ~ 1. In contrast, with probability one Cis+ r ~ 6xp(— r/£typ) where 
^typ denotes the typical correlation length. It is distinctly different from £ a „ (£typ £,av), having a different critical 
exponent. The localization length is useful because it is proportional to the typical correlation length for quantities 
that are diagonal in the fermion representation.!^ 



The localization length is obtained from integrating equation (64) 

1 D 



X(E) v 

In the commensurate case the following approximation holds for small u 



+ constant. (119) 
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(120) 



Equations ( p2| ) and (|9^) give, for small energy, 

A(^)= ^ ln ^ ^ (i2i) 

I D 111 2D ° u - 

The localization length is infinite only when 5 = and E = 0. 

For both the pure system and the random system, £ tJ ,p ~ A(0)~ 1 ~ A -1 , indicating that i^ yp = 1 and that this 
critical exponent is not modified by the presence of disorder. This result also agrees with the RSRGDT. Balents and 
Fisher studied the same Dirac equation and examined the decay of the average Green function. Hence, they found 
the critical exponent associated with the average correlation length, £ av ~ A~ 2 . 



VII. PROPERTIES OF THE INCOMMENSURATE SOLUTION 



It was shown in section tTV B 2| that the incommensurate solution with 8=^=0 describes the XX chain which has 



no anisotropy. We can use equatit 
or a number of other methods ,E3t 



to find the density of states. Alternatively, using Fokker-Planck equations, 
f it can be shown that the number of states below E is, 



N(E) 



Dpo sinh(7ra) 

7T |J iu (<J)| 2 ' 



(122) 



where u = E/D. The density of states, p, may be found by taking the derivative of N(E). The density of states for 
small u is, 



1 



When 8 = 0, for any u, 



P(E) 



p(E) = p Q . 



(123) 



(124) 



Hence, the incommensurate density of states is always finite. When 8 = the density of states is constant. Thus 
for an XX random chain in a non-zero transverse field there is no Dyson singularity. This agrees with the results of 
Smith.El 

To find the localization length we take the real part of equation (Eij) and integrate over E. For small u, the 
localization length for the incommensurate case is, 



A8h(S) 
Io(S) 



(125) 



The constant of integration must be evaluated by deriving the localization length from other meth.ods.E3 The most 
important property of this result is that unlike for the commensurate case the localization length is always finite. This 
means that the typical correlation length of the corresponding spin model does not diverge when the pure system is at 
criticality. Hence, in a non-zero transverse field the anisotropy phase transition does not occur if there is randomness 
in the transverse field or the isotropic exchange. 



VIII. CONCLUSIONS 



We presented some exact results for the effect of disorder on the quantum critical properties one of the simplest 
models to undergo quantum phase transitions: an anisotropic XY spin chain in a transverse field. By taking the 
continuum limit of the corresponding non-interacting fermion model we were able to map various cases of the model 
onto a a Dirac equation with a random mass. This mapping has the distinct advantage that a number of different 
techniques can then be used to obtain exact analytic results for the density of states and the localization length. In 
the presence of disorder Ising transition of the model is in the same universality class as the random transverse field 
Ising model. If there is only randomness in the anisotropy then the anisotropy transition is also in this universality 
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class. However, if there is randomness in the isotropic part of the exchange or in the transverse field then in a non-zero 
transverse field the anisotropy transition is destroyed by the disorder. By examining the energy dependence of the 
density of states we showed that the dynamical critical exponent, show the existence of a Griffiths' phase near the 
transition, and show that the ground state energy has an essential singularity at the transition. The results obtained 
for the typical correlation length, the dynamical critical exponent, the finite-size scaling of the energy gap, and for 
the temperature dependence of the specific heat near the Ising transition agree with the results of the RSRGDT and 
numerical work. Since our result is explicitly exact this agreement is consistent with Fisher's claim that the RSRGDT 
gives exact results for critical behaviour. The real challenge is whether the mapping to the fermion model used here 
can be used to obtain results for distribution functions and spin correlation functions. Recently soi^Jias^bjeen done 
on distribution functions associated with the zero energy eigenstates of the random Dirac equationoo oLj 
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TABLE I. Experimental realizations of random spin chains. These chains all involve antiferromagnetic exchange except for 
SraCuPti-^Ir^Oe which involves both ferromagnetic and antiferromagnetic exchange. 



Material 



Spin per site 



Reference 



| 

2 

alternating 
l 

2 
1 
1 

2 



Qumolinium(TCNQ)2 
SraCuPti-^Ir^Oe 
MnTPP-TCNE(solvent) 
MgTiOB0 3 (warwickite) 
MgVOB0 3 (warwickite) 
Cu(3-methylpyridine)2Cl2 



and 2 



TABLE II. Critical exponents for the transition in the transverse field Ising chain, without and with disorder. A is a measure 
of the deviation from the critical point. The exponents for the random case were calculated by Fisherliij Some exponents are 
expressed in terms of the golden mean, 4> = ^(1 + H is an external field in the x direction, i.e., the same direction as the 
magnetization. The exponents a and 7 are not defined (n.d.) in the random model due to the presence of the Griffiths phase. 
With disorder S = 00 because (of) ~ (ln(l/|#|))- /5 . 



Exponent Definition No disorder with disorder 

a e ~ A 2 ~ a 0+(log) nxT 

P K>~A" 1/8 2-c/> 

7 X„~A^ 7/4 n.d. 

5 «) ~ H 1 ^ (A = 0) 15 00 

v £ av ~ A - " 1 2 



(cjrCTo) ~ J" 1-17 (A = 0) 5/4 0-1 

T ~ ^ 1 OO 



J x + J>' 





FM x V ^ 


\ 

\ 

\ 





h 


1, 






/ 

/ 




FM y f . 











FIG. 1. Phase diagram for the anisotropic XY spin chain in a transverse field at zero temperature and in the absence of 
disorder. The heavy lines represent second order phase transitions. The horizontal line will be referred to as the anisotropic 
transition and the vertical line as the Ising transition. PM denotes a paramagnetic phase and FM^ denotes an Ising ferromagnet 
with magnetization in the x direction. To the right of the dashed line the energy gap in the excitation spectrum always occurs at 
the Brillouin zone boundary (compare Fig. |^) . To the left of the dashed line gap occurs at a wave vector that is incommensurate 
with the lattice. 
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FIG. 2. Typical dispersion relations for the excitation spectrum of the Hamiltonian (|lj) in the absence of disorder. The two 
cases shown correspond to when the quantity a, defined in equation (p3|), is (a) larger than one and (b) less than one. Note 
that for (a) the energy gap always occurs at the Brillouin zone boundary whereas for (b) it occurs at a wave vector that is 
incommensurate with the reciprocal lattice vectors. The cases (a) and (b) occur in regions of the phase diagram to the right 
and left of the dashed line in Fig. Fl (The commensurate case also occurs on the vertical line h = 0: then fco = tt/2). 
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FIG. 3. Phase diagram of the Ising transition in a random transverse field. The horizontal axis is a measure of the deviation 
from criticality in the non-random model. The vertical axis is the amount of disorder. The four phases are: ferromagnet(FM), 
weakly ordered ferromagnet (WO-FM), weakly ordered paramagnet (WO-PM), and paramagnet (PM). The weakly ordered 
phases are Griffiths phases in which the linear susceptibility diverges but there is only short-range order. Note that the dashed 
line does not represent a true phase transition and that higher order susceptibilities will diverge in larger regions of the phase 
diagram. 




FIG. 4. Energy dependence of the disorder-averaged density of states for the commensurate case for various values of 
the dimensionless parameter 8 (see Eq. (0)), which is a measure of the deviation from criticality. The density of states is 
symmetrical about the Fermi energy (E — 0) and diverges at (E = 0) when 8 < \. This parameter range corresponds to a 
Griffiths phase. Note that only far from criticality (8 3> 1) is there effectively a gap in the system. This contrasts with the 
disorder-free case, for which there is always a gap except at the critical point. 
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